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■ ' We demonstrate the danger of a naive identification of kinetic coefficients in the free-electron 

^ ^ , response expansion in powers of an external electric field and coefficients in the expansion of the 
corresponding equilibrium quantity in powers of the fluctuating intrinsic electric field. A particular 

I ' example of the direct current generation by an AC electric field in a mesoscopic system of non- 

^ , interacting electrons vs. the Fock contribution to the equilibrium persistent current in a system 

' of interacting elect rons is considered. We comment on a recent paper by Kopietz and Volker, 

Tij- I |3ond-mat/9711226| . 



The problem of an adequate theoretical description of electron-electron interaction and electron relaxation processes 
in disordered mesoscopic conductors has become a very hot topic nowadays The complexity of the problem 

^ stimulates the search for approaches which incorporate general relationships between kinetic and equilibrium properties 

of the system. One may believe that the knowledge of the nonlinear response to an external field supplemented by 
"Y"; the field correlation functions obtained, e.g., with the use of the fluctuation-dissipation theorem (FTD) ||^ is helpful 
for solving the problem. 

^ The purpose of this Letter is to demonstrate how dangerous might be a naive way of applying FDT-like relations 

' for description of systems of interacting electrons, when the coefficients of expansion of an equilibrium quantity in 
powers of the fluctuating intrinsic electric field are identified with the kinetic coefficients in the free-electron response 
Q ' expansion in powers of an external electric field. 
Q , The reason is a deep physical difference between the dissipative nature of the causal response to an external field 
and the "non-causal" nature of all the interaction processes in the equilibrium. 
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As an actual example we consider the nonlinear generation of the direct current Ijj^ of non-interacting electrons 
in a mesoscopic ring of the circumference L pierced by a constant magnetic flux (f> and an AC flux (f)i^ that produces 
• an electric field £{uj) = [iuj / Lc)4>i^ . This kinetic problem has been studied earlier in Ref. Q and an expression for 
the nonlinear conductance ^^^^(a;, —lo) defined by 

has been obtained. A peculiar feature of the result of Refs. Q is that the nonlinear conductance cr^^,^(w, —lo) does 
not vanish exponentially in the limit of high frequencies uj ^ E^, where Ec = D/L^ is the Thouless energy and D is 
the diffusion coefficient. This is in a striking contrast to a similar mesoscopic effect of the Aharonov-Bohm oscillations 
of the linear conductance 

Recently there was an attempt by Kopietz and Volker |^ to relate the nonlinear conductance cr^i^(w, —u) with the 
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coefficient aeq {q_, ia;„) in the expression for the Fock (exchange) contribution to the equilibrium persistent current in 
the system of interacting electrons: 

AF) _ ^V^,(2)/„ . , ((q-Eq(ia;„))(q-E_q(-za;„))) 

^PC - 1^ Z^'^eq XH^^^n) |— p • l^J 



q.n 



Here ujn are bosonic Matsubara frequencies w„ = 27rTn, and {E'^{iujn) E_^{-~iujn)) = qaqp D{q^,iuJn) is proportional 
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to a disorder-averaged Green's function of longitudinal photons ||^ D(q, zcj„) = ATr/q'^e{q, i|w„|) in a media with the 
dielectric function e(q, w). 

One can transform Eq.(0) to an integral over real frequencies by representing the sum in terms of a contour integral. 



AF) _ St '^j{F)( 

lp(j - ipc + 27r ^ 
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where Ip(j is a "static" contribution given by the a;„ = term of Eq. and 



/lfj(q,c.) = coth (^) [za(2)(q,^)'z'^(q,c^)] ■ (4) 

Here (Teq{c[,uj) — [aeq {—c[, —uj)]* and D{c[,uj) = [£>(— q, — a>)]* are analytical continuations of (Teq {q, iojn) and 
£'(q, ia;„) in the upper half-plane of to. 

Eq.(^ can be further simplified in the Ohmic regime characterized by a constant conductivity cto, where e{q,uj) = 
4:T:aa/{—iuj). In this case _D(q, cj) is purely imaginary and is related with the spectral density of the fluctuating 
intrinsic electric filed (|Eq,^p) by the FDT |]: 

(|Eq,^|2)=zcoth(^) q'D{ci,u;)=coth[^) ^. (5) 

Thus we arrive at an expression analogous to Eq.(jl|): 

jjfj(q,c.)=5?a(2;(q,c.)(|Eq,j2). (6) 

It is tempting to interpret Eq.(^ as a "nonlinear response" to the fluctuating intrinsic fleld Eq^(^ in the system of 
interacting electrons, by an analogy with the response Eq.(|l|) to an external field £q^^ applied to a free electron system. 

Then one may believe |^ that there is a simple relationship between cr^^)j('^, — i^) in Eq.(0) and (Teq^(q, lo) in Eq.(^: 

= limq-O (q; CO). (7) 



However, this leads to a dramatic discrepancy with Ref. 



3J3], since aiq\q]uj) appears to be exponentially small 



at w £'c, while the cr^j^^ni^^, — w) decreases only as u; ^ l^:^!- 

Below we will show that the relationship is principally wrong. This means that one cannot obtain a correct 
Fock contribution to the equilibrium persistent current by substituting Kffeg (q; w) for the free-electron nonlinear 
conductance cr^^^(q, —uj) (at a finite momentum q) into Eq.(|6|). 

To this end we use the expression ||^ for deq^ (q, iu;„) as the derivative of the free-electron polarization operator 
no(q, iitJ„) over the magnetic flux (j): 

cr(2)(q, w„) = ^^Ho(q,icj„), (8) 

where the overline denotes the disorder average and 

Ho(q, iw„) = r ^ Gk+q.k'+qlie; + iw„)Gk',k(iei) = T^Tr^Giiei + iuj^) G{iei). (9) 

k,k',i I 

Here G are Matsubara Green's functions of non-interacting electrons in the presence of disorder; ei = {21 + 1)ttT, and 
the operation Trk is deflned to shorten the following expressions. 

The analytical continuation in Eq.(^) can be performed in a general form: 

Ho(q,c.)=rrk 1^1^ [[G% + c.)G%)-G^(6)G^(e-c.)]tanh(^) + (10) 
+ 2G^(e + tj)G^(e) [^^(e) - ripie + lj)]] , 

where np{e) = [exp (e/T) + 1]^^ and the retarded (i?) and advanced {A) electron Green's functions are the analytical 
continuations of the Matsubara Green's function G(e) into the upper (Ime > 0) and lower (Ime < 0) half-planes, 
respectively. 
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The flux derivative corresponds to breaking of one Green's function into two by the current vertex j: 



(11) 



It is important that neither the analytical properties nor frequency is changed in the Green's functions involved. 
Thus, the function in the r.h.s. of Eq.(|^ takes the form: 



3?a(2)(q,o,) ocTrk 



de 
Atti 
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+ G^(e + c.)G^(e)jG^(e) - G^{e)G^{e ~ uj)jG% - u;)] tanh (^) + 

+ 2[G^(e + w)jG«(e + a;)G^(e) + G^(e + c^)G^(e)jG^(e)] [npie] - npit + w)]] + {uj ^ -uj). 

Eq.(|l2|) can be also obtained by the direct analytical continuation of the expression K{iLUi, icum) given by the triangle 
diagram in Fig. la.: 



K{iu}i,iujra) = Tri^T^G{ien)G{ien + iuji)jG{ien - iuj„i). 
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Fig.la. 
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Fig.lb. 



FIG. 1. a). A generic traingle diagram corresponding to Eq.(p^. Bold lines denote electron Green's functions with Matsubare 
frequencies indicated; b).The Fock contribution to an equilibrium persistent current. The field vertices are connected with a 
photon Green's function. 



The equilibrium Matsubara diagrammatic technique (see Fig.lb.) requires iuii — —iujm, since the corresponding 
vertices are connected by the photon Green's function D(q, iuji). That is why in order to obtain Eq.(p^) one should 
first set iuJi = —iuim in Eq.(p^ and then perform an analytical continuation of K{i(jji) = K{iuji, —iui) into the upper 
half- plane of oj. One can check that: 

3fia(2) (q, uj) (X K^iiu) + K^{-uj) + c.c. (14) 

The nonlinear conductivity (T*-^^(a;, —lo) can also be obtained from a generic triangle diagram Fig.la. However, the 
proper procedure is completely different. It is determined by the causal nature of the system response to an external 
field. For a generic quadratic response to a time-dependent external field E{t): 

/oo /"OO 
dn / dT2<y{Tl,T2)E{t-n)E{t-T2). (15) 
-co <y —oo 

causality requires that (t(ti,T2) oc 9{ti)9{t2), where 9{x) is a step function. This means that the generic nonlinear 
conductivity a{uji,uj2) is an analytical function in the upper half-plane of both lo\ and 0^2 • 

This is the key point in obtaining the nonlinear response from the Matsubara diagrammatic technique j^. One 
should first find an analytical continuation X'^^''^' (wi, ^2) of a generic triangle diagram K{iuji,iuJrn) into the upper 
half-plane of both frequencies, and then set uji — —0^2 = a; to obtain a DC nonlinear response: 
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(2) / \ 
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By performing the analytical continuation in Eq.(|T^) we arrive at: 

de_ 

2[G^(e)G^(e + Lu)jG^ie + uj) - G^(e)G^(e + Lo)jG^{e + io)][nF{e) - npie + lo)] + (uj ^ -lo). 



[G^(6 - Lo)G^{e)jG^{e) - G^{e - Lo)G^ie)jG^ie)] tanh (^^) + 



(17) 



Now it is explicitly seen that the R — A structure of Eqs.(|l2|),(|r^) is different even before disorder-averaging is done. 
The quantity ai'^^ (q, uj) in Eq. (|lj) has a structure of a derivative of the linear AC polarizability with respect to the 
DC flux. Therefore the there is no change of analytical properties in passing the current vertex. On the contrary, the 
nonlinear conductance cr^j^ (w, —lo) is a derivative of the AC conductance with respect to the AC flux. Correspondingly, 
R is changed to A in passing the current vertex. This proves invalidity of the basic assumption @) of Kopictz and 
Volker Ref. The quantities Eqs.([l^),(|l7|) have different physical meaning and the difference in their analytical 
structure is crucial for the functional dependence of their disorder-average values at high frequencies. 

To this end, we note that only such terms in Eqs.([l^),(p7|) that contain both retarded (R) and advanced (A) electron 
Green's functions make a significant contribution to the disorder-averaged values. This is because of the diffusion 
propagators (diffusons and cooperons) P(q, w) oc TrqG^{e -\- Lo)G^{e) which are singular at small q and lo. It implies 
that only last lines in Eqs.( p^,(|l7| ) are relevant. Because R is changed to A and frequency is unchanged in passing 
the DC current vertex in Eq. (|17D, it is possible to build a diffusion propagator (cooperon) of zero frequency in the 
kinetic case. There is no such a possibility in the equilibrium case, since the R-A structure is unchanged in passing 
the DC current vertex. The relevant two-cooperon diagrams are shown in Fig. 2. 
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Fig.2a. 



Fig.2b. 



FIG. 2. a). The leading contribution to the disorder-average kinetic coefficient (7^^)^{lo , —lo); b). The leading contribution 
to the disorder-average value of creq'(0, tj). The bold lines denote electron Green's functions; the shadowed area between two 
dotted lines denotes the fan series (cooperon) P(q, lo) with the frequency lo equal to the difference of frequencies of electron 
Green's functions which are connected by the dotted lines. 



It turns out that both terms in the last lines of Eqs.(p^),(pT[) lead to identical contributions to cr^^'j(tj, — cj) and 

SftCTeg (OjO;), respectively. In the case of nonlinear response where one of the cooperons is of zero frequency, the 
summation over the cooperon momenta is mainly contributed by small momenta \k\ ~ 1/L <C l/i^^ = \fuojD. 
Then one can set fc = in the second cooperon and obtain a power-law cj-dependence. In the equilibrium case 
both cooperons are of the same frequency w, since they appear due to differentiation with respect to the DC flux of 
a single-cooperon diagram for the Aharonov-Bohm effect in the linear AC polarizability. As a result, the diagram 
turns out to be exponentially small in full agreement with the exponentially small Aharonov-Bohm oscillations of the 
polarizability at high frequencies. 

In conclusion, we have shown that the analytical structure of the expression for the free-electron nonlinear con- 

(2) (2) 

ductance cr^j^(q, lj, — cj) is different from that of the coefficient aeq {<1',lo) that couples the Fock contribution to the 
equilibrium persistent current and the spectral density of the electric field fluctuations in the system of interacting 
electrons. The difference is caused by the causality requirement that determines the nonlinear response to an exter- 
nal field but is irrelevant for the interaction contribution to equilibrium quantities. As a result, at high frequencies 
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u) ^ Ec the free-electron nonlinear conductance al.J^{uj , —uj) oc l/w^, while Ueg (q — ^ 0;ci;) is exponentially small. 
This analysis warns against a possible mistake when one applies methods developed for describing kinetic phenomena 
to equilibrium systems and vice versa. 
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